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1. Introduction. The standard theory of the tearing-mode evolution identifies three stages. 
The first is the hnear stage described by the Furth-Killeen-Rosenbluth (FKR) theory pp. 
During this stage, the island width W grows exponentially in time until it reaches the width 
of the resistive dissipation layer, £^ oc r/^/^A'^/^, where rf is the resistivity and A' is the 
instability parameter [see Eq. Q]. Once W ~ nonlinear terms are sufficiently large to 
replace inertia as the force opposing the infiow pattern. A slow down of the growth ensues, 
from exponential to linear in time: dW/dt ~ rjA' . This is the second stage of the tearing- 
mode evolution, known as the Rutherford regime |2|. Finally, the third, saturated, stage is 
reached when the island width becomes comparable to the equilibrium shear length [3]. 

In this paper, we find the tearing-mode evolution to be, in fact, a four-stage process: 
the FKR regime, the Rutherford regime, a regime of fast nonlinear island growth that we 
identify as Sweet-Parker (SP) reconnection, and saturation. We carry out a set of numer- 
ical simulations that demonstrate two main points. First, we show that, given sufficiently 
small ?7, the Rutherford regime always exists; larger values of A' require smaller values of r]. 
Rutherford's negligible-inertia assumption is validated and the asymptotically linear depen- 
dence of dW/dt on 7] and A' is confirmed. Second, we find that, at large A', the Rutherford 
regime is followed by a nonlinear stage of fast growth linked to X-point collapse and forma- 
tion of a current sheet. This causes the reconnection to become SP-like. The signature rj^f^ 
scaling of the effective island growth rate is, indeed, found in this nonlinear stage. The SP 
stage culminates in the saturation of the mode, which can, thus, be achieved much faster 
than via Rutherford regime. 

2. The Model. We use the conventional Reduced MHD set of equations |lj in 2D for a 
plasma in the presence of a strong externally imposed magnetic field 3^. 
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— + v^-Vuj = B^-Vj\\, {1] 



^^+v^-Vij = 7]V'ij, (2) 

where the total magnetic field is B = B^e^ + -B^, all gradients are in the (x,?/) plane, 
the in-plane magnetic feld is B_l = x V^/^, the in-plane velocity is = x V0, and 
the parallel components of the vorticity and current are uj = ■ (V x i?) = and 
j|l = ■ (V X B) = V^'?/'. Eqs. are solved in a box with dimensions x Ly. All 

lengths can scaled so that the width of the box is = 27r. 

We impose an initial equilibrium defined by ip^^^ = V^o/ cosh^(x) and (p^'^^ = 0. We 
choose ipQ = 1.3 so that the maximum value of the unperturbed in-plane magnetic field 
X is 1. Time is scaled by the in-plane Alfven time. The equilibrium is 
perturbed with tp^^^ = ipi{x) cos{ky) , where k = mLx/Ly. In our simulations, the initial 
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FIG. 1: Scatter plot of j|| vs. ip. Run parameters 
are A' = 8.15, 77 = 2.8 x 10"''. Data extracted at 
W = 0.9 [of. Fig. 121. At the separatrix, tp ~ 1.17. 



FIG. 2: Island width vs. time for A' = 8.15 and 
several values of 77. 



perturbation has m = 1. The island width is then approximately W = 4a/ ^/'i(0)/'?/'(*^)"(0). 
However, in what follows, W is measured directly from the numerical data. 

With the equilibrium configuration we have chosen, the instability parameter is 0: 



^i(O) [ki{kf-4) 



6A;2-9 



(3) 



where kf = k'^ + 4. The equilibrium is tearing unstable if A; < y/b. We vary the value of A' 
by adjusting the length of the box Ly. 

3. The Rutherford Regime. Rutherford's analysis depends on the assumption of negligi- 
ble inertia, which reduces the vorticity equation to -B^- Vj||=0. This implies that j\\ = j\\{ip) 
everywhere except at the separatrix (the in-plane magnetic field vanishes at the X-point). 
Fig- E proves the validity of this assumption. At the separatrix, ip ~ 1.17. Larger values of 
ip correspond to the interior the island. The variation of ip in that region is relatively small, 
supporting the "constant-?/' approximation" used by Rutherford. 

Fig. El shows the time evolution of the island width W at constant A' = 8.15 and varying 
resistivity. After the exponential growth stage, a distinct period of linear in time growth is 
manifest in all curves. Figs. 01 and |^ show the dependence of dW/dt during the linear in time 
period on A' and t], respectively, demonstrating in both cases the linear relation predicted 
by Rutherford. At fixed finite rj, the Rutherford scaling breaks down at large A'. However, 
for a given A', it is recovered asymptotically at sufficiently small 77. 

Although Rutherford-like island growth has been observed in earlier numerical work, no 
parameter scan showing the linear scaling of dW/ dt with r] and A' has previously been per- 
formed. White et al. ^ verified the linear in time growth of the island in their numerical 
simulations of the m = 2 mode performed in tokamak geometry, a result later confirmed by 
Park et al. [7 . Biskamp jH] demonstrated the Rutherford behaviour in a numerical experi- 
ment done in slab geometry and with a current-dependent resistivity. Biskamp's simulations 
had a relatively small value A' = 3. Recently, Jemella et al. [0] carried out a A' parameter 
scan with A' e [0.92,20.93] (and constant rj). Their results cast doubt upon the validity 
of Rutherford's analysis by failing to produce the linear in time behaviour: except for the 
smallest values of A', the island growth was exponential at all times. They argued that, 
instead of Rutherford's X-point configuration, a Sweet-Parker (SP) current sheet is formed 
and, consequently, SP reconnection replaces the Rutherford regime. 
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FIG. 3: dW/dt vs. A' at constant = 2.8 x lO""'. FIG. 4: dW/dt vs. 77 at constant A' = 4.97, 8.15. 



We think that the absence of the Rutherford regime in Jemella et al. simulations is, in 
fact, due to the particular equilibrium configuration that they used: B^'' = cos(x). We 
have also performed simulations with such a configuration and found the Rutherford stage 
absent, validating their results. The reason for this apparent sensitivity to the equilibrium 
profile is that a cosine configuration gives rise to two islands in the simulation box. In the 
nonlinear stage, circulating flows between the two X-points are established, which impedes 
the formation of well-separated nonlinear slow-down flows derived by Rutherford. We believe 
that, in order for the Rutherford regime to be recovered in the cosine conflguration, even 
smaller values of rj than used by Jemella et al. (or by us) are necessary. Biskamp who used 
By^^ = tanh(x), did not have this problem and, accordingly, was able to see the Rutherford 
regime. 
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FIG. 5: Contour plots of the current ji\{x,y) at before (t=352, left panel) and after (t=385, right panel) 
for a run with A' = 12.2, rj = 2.8 x 10^"* the nonlinear speed up (the same run as in Fig.EJ. 

4. The Fast Nonlinear Stage. In simulations with large A', we flnd that the X-point 
conflguration maintained during the Rutherford regime eventually collapses and a current 
sheet is formed (Fig. Ej). A dramatic speed up of the island growth ensues (Fig. Fig. [3 
shows that the peak effective growth rate 7 = dlnW/dt ~ 1/t during this fast nonlinear 
stage scales as t]^^'^, thus supporting the expectation that the reconnection in this regime 
is of the SP kind. This behaviour can be qualitatively understood by recalling that the 
X-point conflguration assumed in Rutherford's analysis is, in fact, unstable [lOj. As time 
goes on, the effective island growth rate in the Rutherford regime decreases (7 ~ 1/t) and 
eventually becomes so low that the X-point conflguration can no longer be sustained over 
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FIG. 6: Effective growth rate 7 and island width 
W vs. time, for a run with A' = 12.2 and 7] = 
2.8 X 10"*. 



FIG. 7: Dependence of the peak effective growth 
rate 7max in the fast nonhnear regime on 77^/^ at 
constant A' = 17.3. 



the time scale ~ 7"^. X-point collapse leads to the formation of the current sheet. 

Note that, for large A', if 77 is not sufficiently small for the Rutherford stage to occur, 
the SP stage nearly immediately follows the FKR stage, as is the case in Jemella et al. 
simulations. In the opposite case of small A' and small ?], the saturation can be reached 
directly from the Rutherford stage, with the SP stage never having time to materialise. 

We leave a more detailed theoretical description of the nonlinear speed-up effect to a 
forthcoming paper. 
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